32-goyal

Elio Campitelli

Exploración de los datos de Goyal et al. (2021)Goyal, Rishav, Martin Jucker, Alex Sen Gupta, Harry H. Hendon, and Matthew H. England. 2021. “Zonal Wave 3 Pattern in the Southern Hemisphere Generated by Tropical Convection.” Nature Geoscience, August, 1–7. https://doi.org/10.1038/s41561-021-00811-3..

En el paper hacen simulaciones globales con distintos continentes y ven qué le pasa a la onda 3. En el repositorio con los datos están disponibles las simulaciones:

Faltan la simulación control, la que tiene sólo la parte subtropical de Sudamétrica, la que tiene la parte subtropical de todos los continente, y una que tiene sólo la Antártida.

0.1 Ondas zonales

La Figura 1 reproduce la Figura 2 del Goyal et al. (2021)Goyal, Rishav, Martin Jucker, Alex Sen Gupta, Harry H. Hendon, and Matthew H. England. 2021. “Zonal Wave 3 Pattern in the Southern Hemisphere Generated by Tropical Convection.” Nature Geoscience, August, 1–7. https://doi.org/10.1038/s41561-021-00811-3.. Muestra las ondas 3 anuales medias para los 100 años de cada simulación junto con la onda 3 media total en rojo.

Figure 1: Onda 3 anual media de cada simulación y para ERA5.

Onda 3 anual media de cada simulación y para ERA5.

Se reproduce bien esto que los continentes están bloqueando la fase de la onda 3, pero que la amplitud media de la onda 3 es grande en todas. En particular, la simulación con Australia tropical tiene una onda 3 enorme, pero en cualquier lado. Pero aún el aquaplanet tiene onda 3 no despreciable.

La Figura 2 siguiente muestra la distribución de amplitud de las ondas 1 a 4. Se ve que todas las simulaciones con algún continente tienen una distribución virtualmente idéntica de ondas 2 y 3; la única simulación con una distribución distinta es la aquaplanet, que tiene amplitud algo menor. La onda 1 es más similar en todas las simulaciones; hay algunas diferencias, pero son pequeñas en comparación con la diferencia con los datos observados. La onda 4 también es casi igual en todas las simulaciones.

Figure 2: Distribución de la amplitud mensual de la ondas 1 a 3 en 55ºS para las distintas simulaciones. En negro, lo mismo para ERA5.

Distribución de la amplitud mensual de la ondas 1 a 3 en 55ºS para las distintas simulaciones. En negro, lo mismo para ERA5.

La distribución meridional de la amplitud media de las ondas también es prácticamente idéntica en todas las simulaciones con algún continente, como se ve en la Figura 3. De nuevo hay algunos cambios en la onda 1, que en la simulación con trópicos realistas alcanza una amplitud comparable a la observada pero más al sur.

Figure 3: Amplitud media de las ondas zonales por latitud.

Amplitud media de las ondas zonales por latitud.

En resumen, para la onda 3, la simulación aquaplanet tiene una amplitud ligeramente menor, pero no lo suficiente para explicar el cambio dramático en la amplitud de la onda 3 del campo medio. La diferencia entre simulaciones se encuentra en la estacionariedad de la onda.

La Figura 4 cuantifica la estacionariedad de la onda 3 como la amplitud de la onda promedio sobre la amplitud promedio de las ondas mensuales.

Figure 4: Estacionariedad de la onda 3 para cada experimento y ERA5 (línea negra).

Estacionariedad de la onda 3 para cada experimento y ERA5 (línea negra).

Tanto el Aquaplanet como la simulación con Australia tropical tienen una estacionariedad casi nula. De hecho, esta segundaa simulación tiene aún menos estacionariedad en latitudes medias que la primera. Luego, el resto de las simulaciones con algún continente tienen una estacionariedad más alta. Aunque todas las simulaciones tienen menos estacionariedad que la observada en ERA5.

La Figura 5 resume todo esto, mostrando la amplitud media de las ondas zonales en el eje horizontal y la estacionariedad en el eje vertical. Las líneas negras marcan los valores de amplitud de las onda media correspondiente.

Figure 5: Amplitud media de las ondas y estacionariedad de cada experimento y de ERA5 (punto negro) para los números de onda 1 a 4. Línes de contorno marcan la amplitud de la onda media

Amplitud media de las ondas y estacionariedad de cada experimento y de ERA5 (punto negro) para los números de onda 1 a 4. Línes de contorno marcan la amplitud de la onda media

En esta figura se resume también el comportamiento de las distintas ondas zonales. Para la onda 1, los distintos experimentos varían tanto en amplitud media de las ondas como en la amplitud de la onda media. Y aunque ninguno de los experimentos se acercan a las amplitudes observadas, la estacionariedad de la simulación con trópicos realistas (~0.7) es bastante alta y cercana a la observada.

En la onda 2 la cosa es distinta. Muchas de las simulaciones consiguen una amplitud media y una estacionariedad más o menos realista o incluso mayor a la observada.

Finalmente, para la onda 4 todas las simulaciones tienen una onda 4 con amplitud media y estacionariedad comparable a aquaplanet. La excepción es la de trópicos realista, que tiene mucha más estacionariead.

Esto es consistente con lo que vi en la tesis de licenciatura. Las condiciones de superficie tienen el rol principalmente de cambiar la estacionariedad de la onda, no tanto la amplitud. Las ondas se excitan aún en ausencia de continentes. Poner un poquito de continente (Australia tropical) excita aún más, pero no es suficiente para bloquear la fase significativamente.

0.2 ¿Sirve fourier?

Una pregunta es qué tan bien funciona el paradigma de Fourier en estas simulaciones.

Figure 6: Envolvente media de las ondas zonales de cada experimento

Envolvente media de las ondas zonales de cada experimento

Podemos preguntarnos qué tan “concentradas” están las ondas en cada experimento. Es decir, si suelen ser ondas planetarias o están más localizadas. La Figura 6 de abajo muestra la envolvente media de las ondas zonales para cada experimento. En todos los experimentos, y en las observaciones, hay un anillo con amplitud de las ondas máxima. En Aquaplanet, este anillo está en latitudes más bajas que en el resto de los experimentos y en las observaciones y es perfectamente zonalmente simétrico.

Figure 7: Envolvente media de las ondas zonales de cada experimento. En etiqueta, el desvío estándar sobre el valor medio en 55ºS.

Envolvente media de las ondas zonales de cada experimento. En etiqueta, el desvío estándar sobre el valor medio en 55ºS.

Este cambio en la amplitud de la envolvente se coincide con el cambio en la latitud donde es la amplitud máxima de la onda 1 (Fig. 3), y muy probablemente tenga que ver únicamente con esa onda. Para eliminar ese efecto, la Figura 7 muestra la envolvente media de las ondas con número de onda mayor a 1. En texto blanco, se muestra el desvío estándar sobre el valor medio de la envolvente en 55ºS como una medida del nivel de asimetría.

En las simulaciones con Sudamérica y en la de África tropical, hay un máximo en la amplitud de las ondas al sudeste del forzante. La asimetría al rededor de 0.05 es mayor que la de aquaplanet, pero mucho menor que en las observaciones, las cuales tienen un máximo bien definido a lo largo del sur del Pacífico e Índico. La simulación con trópicos realistas también tiene una asimetría similar a las simulaciones.

Esto muestra que en las simulaciones tienen mucha menos preferencia en la localización de las ondas en un lugar en particular que las observaciones. Lo que no muestra es si los campos mensuales son más simétricos o no. Para eso la Figura 8 muestra el desvío estándar de la envolvente de cada campo mensual. Representa una medida de cuán localizadas están las ondas en promedio sin considerar la localización de las mismas.

Todas las simulaciones muestran una localización similar y comparable a las observaciones. Incluso la simulación de Aquaplanet alcanza valores similares, aunque desplazados hacia el ecuador.

Figure 8: Desvío estándard promedio de la envolvente de las ondas zonales con número de onda mayor a 1.

Desvío estándard promedio de la envolvente de las ondas zonales con número de onda mayor a 1.

En otras palabras, lo que varía en cada simulación no es la localización de las ondas zonales sino el sesgo de esta localización en un lugar particular. En el caso de las observaciones, al sur de Australia y del Atlántico, en el caso de la simulación con Sudamérica completa, al sur del Atlántico.

Este análisis es sobre la amplitud de las ondas en general. ¿Pero qué pasa con la onda 3 en particular?

Se puede calcular la proyección de las anomalías zonales de altura geopotencial estandarizadas sobre la onda 3 asociada ese campo en gajos de 60º de longitud de ancho. Esto da una idea de qué regiones están alineadas y “aportan” más con la onda 3.

La Figura 9 muestra la proyección media. Aquaplanet tiene un campo casi constante, indicando que no hay ninguna región preferencial donde la altura geopotential covaríe con la onda 3. El hecho de que no haya una banda más o menos homogénea de latitudes con valores altos indica que la onda 3 en esta simulación está bastante localizada; es decir que no es una onda hemisférica.

Las otras simulaciones tienen una banda de valores altos entre 40ºS y 60ºS que también tiene máximos locales, indicando menor localización. Lo que es interesante es que en las regiones tropicales hay valores altos sobre las áreas con continentes en cada simulación, lo cual es consistente con que ésta son las fuente de las ondas.

La estructura de esta variable en ERA5 es bastante distinta. Si bien hay valores relativamente altos en latitudes medias, éstos no se distribuyen homogéneamente en una latitud sino que cambian de magnitud según la longitud. En particular, en la región del Índico el máximo se encuentra más al norte que en la región del Pacífico y existe un mínimo local importante al norte del Mar de Amundsen. En los trópicos hay máximos bien marcados sobre el Pacífico ecuatorial.

Proyección media de altura geopotential en la dirección de la onda 3 tomando gajos de 60º de ancho.

Figure 9: Proyección media de altura geopotential en la dirección de la onda 3 tomando gajos de 60º de ancho.

Lo siguiente es evaluar la misma proyección que antes pero sobre la onda 3 climatológica. La Figura 10 muestra la proyección media de la altura geopotencial en la dirección de la onda climatológica 3 de cada mes.

Igual que Figura 9 pero en la dirección de la onda 3 climatológica.

Figure 10: Igual que Figura 9 pero en la dirección de la onda 3 climatológica.

La simulación de Aquaplanet que tiene casi nula estacionariedad, tiene valores casi nulos, lo cual es esperable. Como las ondas zonales 3 no tienen una región preferencial, en promedio, la proyección de la altura geopotential en la dirección de la onda climatológica es cero.

Las otras simulaciones tienen regiones con valores positivos. En las simulaciones con Sudamérica, hay un máximo sobre Sudamérica tropical y a sotavento en latitudes medias. La forma de estos máximos es consistente con el trayecto esperado de la onda 3. En la simulación con África hay valores positivos sobre África que parecen extenderse a latitudes más bajas, pero dicha extensión es mucho menor que en las simulaciones con América. En la simulación con Australia, nuevamente hay valores altos sobre la fuente, pero en este caso no hay ninguna propagación fuera de la región. En la simulación con trópicos realistas, hay dos regiones tropicales con valores positivos: Pacífico central y áfrica.

0.3 cEOFs

Como lo que interesan son las onda zonales, voy a calcular los EOFs complejos. Para ayudar en la interpretación de lo que ya sé, voy a usar sólo Septiembre-Octubre-Noviembre.

Figure 11: Parte real de los primeros 3 EOFs complejos de los experimentos.

Parte real de los primeros 3 EOFs complejos de los experimentos.

La Figura 11 muestra la parte real de los primeros 4 cEOFs de las simulaciones. El primer cEOF es una onda 1 en todas las simulaciones, la cual estaría uniendo los primeros 2 EOFs (Fig. ??). Los siguientes cEOF varían más en estructura entre los distintos experimentos.

Salvo aquaplanet y “Tropical South America,” en el resto de las simulaciones los cEOFs 2 y 3 parecerían tener su amplitud localizada en hemisferios opuestos.

Salvo la simulación de Sudamérica tropical, todos estos cEOF tienen un grado de localización comparable con ERA5.

Tratar de hacer paralelismos entre cEOFs no es fácil porque la tienen diferencias en cuanto a las regiones donde están más activos por construcción de las simulaciones. Para tener una mejor idea de la estructura de ondas zonales independientemente de la fase, la Figura ?? muestra la amplitud de las ondas zonales 1 a 3 de los cEOFs 2 y 3 de cada experimento (salvo aquaplanet).

Figure 12: Amplitud de las ondas zonales en cada latitud para los EOFs complejos.

Amplitud de las ondas zonales en cada latitud para los EOFs complejos.

Es importante ver la distribución de amplitudes en el cEOF2 de ERA5, ya que sería el “benchmark” de lo que debe ser el PSA. Por ejemplo, el cEOF2 de la simulación con “trópicos realistas” tiene una distribución casi idéntica que el benchmark salvo por la onda 1. Cabe notar también que el cEOF3 de esta simulación también es muy parecido salvo la onda 1; de hecho, la varianza explicada de cada uno es casi idéntica (~8%), por lo que claramente estos modos están muy mezclados. Es posible que se trate de dos modos que tratan de refleja una estructura similar pero localizada; el cEOF2 explicaría la parte del Atlántico-Índico mientras que el cEOF3 explicaría la parte del Pacífico.

Algo similar parecería ocurrir en la simulación de Complete South America. En ésta los cEOFs no están mezclados en cuando a varianza expicada, pero ambos son similares, con una mezcla similar de onda 2 y 3. La diferencia principal es que el cEOF2 tiene más amplitud an el hemisferio occidental mientras que el cEOF3 tiene más amplitud en el oriental.

La simulación de Tropical South America es muy distinta a las demás. El cEOF2 es una onda 3 casi pura y el cEOF3 es una onda 2 casi pura.

Una última herramienta para tratar de entender estos modos es calcular las regresiones entre ellos y la altura geopotential completa. Esto brinda una idea más realista de los campos asosicados a cada modo y cada parte. En particular quiero mirar dos cosas:

  1. Si la parte real e imaginaria forman realmente una estructura de ondas.

  2. Si hay estas ondas forman parte de un tren de ondas tipo PSA.

Figure 13: Regresiones de altura geopotential y los cEOFs. Las primeras dos filas son las regresiones con la parte Imaginaria y Real del cEOF2 y las últimas dos, las del cEOF3.

Regresiones de altura geopotential y los cEOFs. Las primeras dos filas son las regresiones con la parte Imaginaria y Real del cEOF2 y las últimas dos, las del cEOF3.

La Figura 13 muestra estas regresiones.

Para Complete SOuth America, el cEOF2 es rarísimo y no veo que tenga una estructura clara de onda. El cEOF3, en cambio, tiene un poco más de pinta. Aunque la region que afecta está lejos de sudamérica.

Para Tropical South America, ambos cEOFs tiene una estructura de onda zonal, pero en ambos casos son ondas zonales bastante puras. Hay una sugerencia de un tren de ondas en la regresión con la parte imaginaria.

En Tropical Africa, el cEOF2 es bastante raro. El cEOF3 sí tiene más pinta de tren de ondas.

En Tropical Australia, el cEOF3 tiene un patrón de tren den ondas bastante bien definido.

El Realistic Tropics tiene un patrón de tren de ondas en cEOF2, aunque con superposición de una especie de SAM con la parte Real.

Table 1: Identificación de cEOFs con tren de onda tipo PSA

exp cEOF
Complete South America cEOF3
Tropical South America cEOF2
Tropical Africa cEOF3
Tropical Australia cEOF3
Realistic Tropics cEOF2
ERA5 cEOF2

Las decisiones de qué cEOFs son los relevantes en cada simulación están en la Tabla 1.